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ABSTRACT 


Quantum physics revolutionized classical disciplines of mechanics, statistical physics, and electrodynamics. One branch of 
scientific knowledge however seems untouched: thermodynamics. Major motivation behind thermodynamics is to develop 
efficient heat engines. Technology has a trend to miniaturize engines, reaching to quantum regimes. Development of quan¬ 
tum heat engines (QHEs) requires emerging field of quantum thermodynamics. Studies of QHEs debate whether quantum 
coherence can be used as a resource. We explore an alternative where it can function as an effective catalyst. We propose a 
QHE which consists of a photon gas inside an optical cavity as the working fluid and quantum coherent atomic clusters as the 
fuel. Utilizing the superradiance, where a cluster can radiate quadratically faster than a single atom, we show that the work 
output becomes proportional to the square of the number of the atoms. In addition to practical value of cranking up QHE, our 
result is a fundamental difference of a quantum fuel from its classical counterpart. 


Superradiance (SR) was originally introduced by Dicke in 1954 as a cooperative emission of light from an ensemble of 
excited two level atoms in a small volume relative to emission wavelength. The atoms radiate in a synchronized (coherent) 
manner, at quadratic rate with the number of atoms. Experimental verification of SR has been accomplished in diverse 
systems. " Due to fast dephasing, typical observations are limited to pulsed or transient regimes, while modern experiments 
achieve SR in steady state as well. Nature herself benefits from the SR processes; in particular, in light harvesting complexes, 
for efficient collection and transfer of solar energy. 0-1 We explore if we can follow the nature’s example and use SR to 
enhance work harvesting capability of a quantum heat engine (QHE). 

Typical QHEs are analogs of classical heat engines, whose working fluid or heat reservoirs are replaced by quantum 
systems. " While the recognition of a quantum system, three level maser, as a heat engine goes back to 1959, many 
modern studies discuss effects of quantum coherence on the performance of QHEs. A particularly intriguing QHE consists 
of an optical cavity pumped by a beam of atoms prepared in coherent superposition states. After repeated passages of the 
atoms one by one, or two by two, the cavity field can be thermalized to an effective temperature which can be higher than 
that of the atom. Using such a quantum coherent thermalization procedure in a Carnot cycle, a QHE with unique properties 
is proposed. “ Photonic Carnot engine could operate with a single thermal reservoir, and could beat the Carnot efficiency. 
Addition of the coherence reservoir and the cost of coherence generation brings these feats within the boundaries of the second 
law of thermodynamics. The work output of the engine is determined by the radiation pressure on the cavity mirrors and it is 
proportional to the cavity field intensity. We propose to enhance the work harvesting capabilities of a photonic QHE using 
SR enhancement of the cavity field intensity. We choose to operate our QHE in Otto cycle instead of Carnot cycle, as it is 
more experimentally feasible. 

We consider clusters of N two level atoms interact with the cavity at regular intervals in the ignition stage of the Otto cycle. 
Steady state SR in an overdamped cluster micromaser ’ is similar to our set up, though we consider weak damping regime. 
The atoms are assumed to be prepared in a low temperature coherent superposition state and the cavity field is assumed to 
be in a thermal state initially. We find that after a number of interactions, the cavity field reaches an equilibrium state, which 
is a coherent thermal state. " The corresponding work output of the QHE can be determined from the steady state photon 
number, which is enhanced by the SR. Accordingly it becomes proportional to N 2 . 

Quantum coherence in the cluster is used here as an effective catalyst, increasing the energy transfer rate from cluster 
to the field. It is necessary to prove that quantum coherence can be completely recycled in order to call it as a true catalyst. 
As there are fundamental and practical challenges against that, " we leave it here as an intriguing open question for future 
contributions if and how quantum coherence can be recovered partially or completely. Our focus here will be to utilize 
quantum coherence to enhance the rate of energy transfer. Coherent cluster acts like a fuel which burns quadratically faster 
than an incoherent one. The time to reach equilibrium increases when the energy kicked into the cavity in rapid bursts. This 
allows for harvesting quadratically more energy from a resource before equilibrium is established. Coherence is partially 
transferred from a set of clusters to the cavity field. The field can be described as a thermal coherent state in equilibrium. 


In contrast to energy, which is harvested as work output, coherence is not consumed and it remains within the field. More 
technically, quantum coherence is determined by the off diagonal elements of the density matrix of the field; and work is 
harvested in a quantum adiabatic stage, where the off diagonal elements are preserved up to a geometric and dynamic phase 
factor. • The diagonal elements, and hence the number of photons, cannot change in accordance with the adiabatic theorem. 
The coherence could in principle be transformed from the cavity field back to the clusters in the exhaust stage. The cluster 
resonance frequency should be altered to match the adiabatic changes in the photon frequency. While one could envision a 
hybrid system where coherence is exchanged between clusters and photons in a completely or partially reversible manner, 
and work is extracted locally from photon subsystem, there are serious obstacles against its realization. Complete recovery of 
initial coherence may not be possible first due to fundamental constraints by the time translational symmetry; “ and second 
due to practical constraints by the rapid dephasing of coherence. In practice one can take precautions against dephasing using 
time dependent control methods, 59 and be content with partial recovery of coherence within the fundamental bounds to reduce 
the production cost of quantum coherent fuel. Instead, coherence can be repeatedly reestablished externally in the clusters 
before every cycle of the engine operation. This is an effective use of catalytic coherence. 


The model and dynamical algorithm 

We consider a dissipative, but high finesse, single mode optical cavity, whose field is interacting with a of cluster of N two- 
level atoms for a time f ml repeatedly at a rate r ^ 1/^int- A beam of clusters passing through the cavity in a time of t ml one at a 
time at regular intervals of 1 jr> f; nt is an equivalent alternative scenario. 

The interaction is described by the Tavis-Cummings Hamiltonian (in units of h = 1) 

77 sys = (Of o'a + (O a S z + g(aS + + a^S~), (1) 


where (Of is the cavity photon frequency, (0 a is the transition frequency of the atoms, and g is the uniform interaction strength. 
The photon annihilation and creation operators obey the boson algebra and denoted by a and a', respectively. The atomic 
cluster is represented by collective spin operators 1/2 £,-of), which obey the St/(2) spin algebra, where 

of 1 and of are the Pauli spin matrices, corresponding to the transition and population inversion operators for the ith atom, 
respectively. The multi photon generalization of the Tavis-Cummings model is considered in quantum Otto cycle from the 
perspective of the interplay between dynamical Stark shift, thermal entanglement and the engine efficiency. 

The dynamics of the state of the system, p sys , during the interaction, is determined by the master equation 
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where K is the cavity dissipation rate. Between the interactions, the cavity field evolves freely by the master equation 
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where P/ = Tr a (p sys ) is the reduced density matrix of the field, and Tr„ is the partial trace over the atomic states. 

Every interaction starts with an initial state of the system in the form p sys = p a ( 0) 0 Pf, where p„(0) is the externally 
prepared state of the cluster, which is the same at the beginning of every interaction; and Pf denote the state of the cavity field 
which is changing from one interaction to the other according to dynamics described by Eqs. (2) and (3). 


Results 

Thermalization 

In a usual set up, the atomic ensemble, first, gets into contact with a hot reservoir at temperature 7},. After a sufficient period of 
time, the state of the individual atoms p‘ a thermalizes, i.e., p‘ a ( 0) = (1/Z«)£„exp(— PE l n )\y/ n y (y/ n \. Here, E„ and \y/ n ) are the 
eigenvalues and the corresponding eigenvectors of the single-atom Hamiltonian 77' = (( o a /2)o ', /l = (1 /7},) (kn = 1) is the 
inverse temperature and Z a = Tr(p/(0)) is the partition function. The state of the cluster, assuming there are no interatomic 
interactions, is given by p a { 0)' = <8>,Pa(0)- The initial state of the cavity field is a thermal state at a temperature T c , i-e., 
Pf(T c ) = (1/Zy)exp(—with Hf = (Ofa*a and Zf = Tr(py(0)). 

We assume that 7), is sufficiently low, such that most of the atoms are in their ground states. Prior to their interaction with 
the cavity field, the atoms are transformed into coherent superposition states by a rotation operation 

p„(0)= J R(C)p a (0)'P(C) t , (4) 
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Figure 1 . Dynamics of the effective temperature and the mean photon number, (a) The effective temperature T e s and 
(b) the corresponding average number of photons (n) as a function of the scaled time cot, respectively, for 100 successive 
interactions of IV = 2 (black-solid), N = 3 (green-dashed) and N = 4 (blue-dotted) atom clusters. Brown-dashed lines 
represent the initial values T c = 0.5 and (n c ) = 0.156. (c) The steady state values of T & g- and (n) as a function of N obtained 
by time averaging over 250 successive injections. All the other parameters are as explained in the text. 
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where R(Q = exp (£S + — C*S~) with £ = (0/ 2)e“ p . In our calculations we set cp = 0 and take 0 = —it/ 2, i.e., the atomic 
states are moved from the pole of the Bloch sphere to the equator. The initial state of the system is given by p sys = p a ( 0) 0 pj. 
Our choice of the initial state is not arbitrary. It is well known that the collective atomic coherent states are closely related to 
the Dicke states and superradiant with the choices of cp = 0 and e> = ±n/2. The initial state of the atomic cluster p a ( 0) can 
be called as a thermal coherent spin state. 

We assume that the field and the atoms are in resonance, i.e., (Of = (O a , which will be used to scale energy, time and 
temperature parameters in the numerical simulations. We take g = 0.19, K = 0.03, 7), = 0.001, T c = 0.5 and rt; nt = 1/6. 
Since g/K > 1, our system is not in the overdamped micromaser regime and we cannot use the corresponding superradiant 
micromaser master equation. ' We solve Eqs. (2) and (3) instead to determine the cavity field density matrix after a number 
of interactions with the clusters of N atoms. Our typical results are shown in Fig. 1. In Fig. 1(a) and 1(b), we present the time 
dependence of the effective temperature T e ff and the mean number of photons in the cavity (n) = Trip/a'a), respectively. The 
initial number of photons in the cavity is calculated to be ( n c } ~ 0.156 at T c = 0.5. The effective temperature is defined by the 
relation (n) = 1 /[exp (1 /7/ff) — 1], It can be interpreted as a temperature only in a steady state which is an approximate thermal 
equilibrium state. If the coherence is large, the probability distribution of the photons in the steady state SR phase becomes 
closer to that of a coherent state, with the increasing N. More rigorously, it can be described by a thermal coherent state. 

It can be written as a thermal state Pf{T) subject to Glauber displacement transformation P/(7’, a) = D(a)Pf(T)I)' (a), with 
the displacement operators D(a), where a is the coherence parameter. 

The results are shown for the cases of clusters with N = 2,3,4 atoms. In each case the time range allows for 100 interactions 
between the clusters and the cavity field. We see that the cavity field takes more time to reach steady state with larger clusters. 
This observation is in fact ensures that catalytic use of coherence is consistent with the energy conservation. Each cluster 
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stores an energy of N(0 a /2, which scales linearly with the number of atoms. If M clusters are needed to bring the cavity field 
to equilibrium, then the energy delivered from clusters to the cavity will be a fraction of MN(O a /2. The change in the energy 
of the cavity field scales quadratically with N. Accordingly, if M would not increase with N , we could in principle transfer 
more energy than the stored amount by simply increasing N to extreme values. Hence, the increase of M with N is consistent 
with the energy conservation. 

Fig. 1(c) shows the steady state values of the effective temperature of the cavity field T e ,g ss and the mean photon number in 
the cavity (n) ss that are calculated via time-averaging over a period of t = 1500 corresponding to 250 successive interactions 
of N = 2,3,4,5,6 atom clusters. Curve fitting yields the relations T e g tSS = T C +0.1N 2 and (n) ss = (n) c + 0.095N 2 . The 
proportionality constant t, ~ 0.1, that appears in front of the N 2 scaling, is due to the single atom micromaser emission 
intensity. In micromasers with random arrival times of pumped atoms, the analytical expression of the emission intensity is 
given by I\ = rg 2 tf m P e / K where P e is the probability of finding the two level atom in its excited state. Our parameters used in 
this formula indeed verifies that incoherent single atom emission is ~ 0.1. On the other hand, we cannot immediately conclude 
that this formula is applicable to the our set up of regular injection of clusters, for which case the analytical verification of 
such a relation is proven to be a difficult problem. Nevertheless, we further test these relations via numerical investigations 
and found that B, 4 g 23 for K = 0.03 and B, 0.008fc °' 7 for g = 0.19. These numerical fits comply with the analytical 
estimation. 

The mean number of photons in a thermal coherent state (TCS) is given by (h)tcs = («)th+ O' | 2 . The effective tempera¬ 
ture description is not essential to comprehend or operate the thermalization. We recognize that the physical temperature of the 
cavity remains the same and the atoms only transfer coherence to the field. Clusters therefore act as a pure coherence reservoir. 
By fitting TCS to the numerically determined density matrix in equilibrium, with fidelity ~ 1, we verify that | a | 2 ~ 0.11V 2 . 
This quadratic enhancement of coherence transfer into the cavity field can be translated into useful work output, by associating 
number of photons with the radiation pressure, in a photonic QHE. 

Superradiant quantum Otto engine 

We consider a four-stroke quantum Otto cycle for our photonic QHE. The working fluid, which is the photons inside the 
cavity, is described by Hamiltonian Hf = (Oa [ a with eigenvalues E n ,(n =0,1,2...). The corresponding eigenstates | n) has 
occupation probabilities P n (T) = exp (—J 3E n ) /Z. 

In the ignition stroke, the photon gas is heated in a quantum isochoric process to temperature 7# and the occupation 
probabilities change to P n (T H ). The eigenvalues remain constant and denoted by E'J. There is only heat intake and no 
work is done. This process is assumed to happen by the superradiant thermalization procedure where the photons interact 
with coherent atomic clusters. For simulations we choose the starting temperature of the cavity field as T c = 0.5, in units 
of (Oh , which is the frequency of the atoms and the field. After the thermalization, the photon gas effective temperature 
becomes Th = 7 d f ss . Both the energy and coherence of the clusters are partially transferred to the photon gas, which can 
be described approximately in thermal state in the case of weak coherence. More generally, we do not need an effective 
temperature description and one can consider larger coherence transfer by thermalization into a thermal coherent state. The 
clusters and the cavity field comes into an equilibrium in terms of coherence, rather than temperature. This vision is in parallel 
with the recent ideas on thermodynamics of quantum coherence " as well as with other generalizations of second law 
and thermodynamical principles. “ If we denote the initial thermal density matrix of the field as Pf(T c ), it evolves into a 
coherent thermal state P/(7’ c , a) in the steady state where | a | 2 ~ 0.1N 2 . 

Second stage is the expansion stroke, which is a quantum adiabatic process, where the work is done by the photon gas 
and the eigenvalues change to E ! n by changing the frequency to ft)/ ; while there is no heat exchange and the occupation 
probabilities remain the same, the physical temperature drops from T c by expansion. The off diagonal elements remains the 
same as well up to a dynamical and a geometrical phase factor. The diagonal elements of the density matrix, and hence the 
photon number, remains the same and the magnitude of the coherence is preserved. Ideally the coherence would not vanish in 
the adiabatic process, but in practice decoherence and dephasing would reduce coherence. The density matrix of the field at 
the end of the expansion can be written as Pf(T\ a’). 

In the subsequent exhaust stroke, the state of the photon gas should be transformed into a thermal one as p (Tl) in another 
quantum isochoric process and the occupation probabilities change to P n (TL), but the eigenvalues remain constant. An amount 
of heat released, but no work is done. Coherence can be transferred to environment by dephasing or in principle to clusters. If 
the same clusters would be used, it is necessary to alter their frequencies to a>i. Due to conservation of the photon number and 
by variation of the atomic frequency, the coherence transfer is energetically allowed, though dephasing and time translation 
symmetry constraints makes it limited. “ Accordingly one could choose more practical methods than using clusters to cool 
the cavity. The thermal density matrix of the field at the end of the exhaust stroke can be written as Pf{Ti). 

Compression stroke is the last stage where the eigenstates change back to Ejf by variation of the photon frequency from 
(Oi to (Oh in another quantum adiabatic process where work is done on the cavity field without any heat exchange; and the 
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Figure 2. Work output (W) of the photonic quantum Otto engine, (a) Work harvested by the photonic quantum Otto 
engine for N = 2 (black-solid), N = 3 (red-dashed) and N = 4 (blue-dotted) atomic cluster pumps, as a function of frequency 
ft)/ (in units of ft)// ), (b) Maximum work output (W max . in units of ft)//) from the engine as a function of N. 


temperature raises back to T c . The density matrix of the cavity field is set back to Pf{T c ). If the same clusters are used in 
cooling stage, their frequency should be changed back to ft)//. Due to dephasing and fundamental constraints on coherence 
transformations, “ the clusters would need to be induced coherence at the beginning of every engine cycle, so that their 
state is set back to p fl (0). The coherence production costs can be reduced by optimizing time dependent control methods to 
reduce dephasing and by using schemes to recover coherence as much as possible. 

In our calculations, we employ the quantum mechanical interpretation of the first law of thermodynamics, where the heat 
absorbed Qi„, the heat released Q out , the net work done W are given by the relations 

Qin = Y, E n[ P n(T H )-Pn(T L )], 

n 

Qou, = Y, E n{ P n( T L)-Pn(T H )], 

n 

W = Qin + (2om = ^ 

n 

where (E ! n ) are the energy levels during the isochoric stages. The efficiency /] is defined by rj = W j Q m . Throughout our 
analysis, we only consider positive work extraction which obeys the relation Q m > —Q 0 ut > 0 in accordance with the second 
law of thermodynamics. 

Using these definitions, it is straightforward to show that the work output is proportional to the difference in the mean 
number of photons W = ri((n) ss — (u)l ), where = («)/, = ( n) c and T] = 1 — (Ol■ These relations are plotted in Fig. 2. The 
work output W is shown in Fig. 2(a) for atomic clusters of N = 2,3,4 atoms. The positive work condition is independent of N 
and given by ft)/, < ft)//. In the simulations, we find the temperature of the photon gas before the compression stroke 7/ for each 
ft)/,, using the relation 7/ = ft)/ 7) which gives 7/ = ft)//2 for T c = 0.5. This relation is a consequence of the requirement to 
retain same occupation probabilities in the quantum adiabatic process. The positive work assumes its maximum at maximum 
efficiency when ft)/ <C 1, accordingly, at low ft)/, photonic quantum Otto engine is able to translate SR enhancement of cavity 
intensity to its work output, as shown in Fig. 2(b). We find that the harvested work is maximum at maximum efficiency and it 
obeys a power law such that W max ~ 0.085A' 2 at T] max ~ 0.99. 

Discussion 

We considered a photonic QHE that undergoes a four-stroke Otto cycle. The working substance is taken as a photon gas in 
an optical cavity. The coherent thermalization of the photon gas has been accomplished via coupling two-level atom clusters, 
acting as quantum coherent fuel, with the cavity at regular time intervals. The cluster atoms have initially been prepared in 
a thermal coherent spin state. The cavity field reaches an equilibrium coherent thermal state. We find that the mean number 
of the cavity photons and the corresponding effective temperature are scaled with the square of the number of the atoms in 
the clusters, due to SR. We argued that the coherence of quantum fuel function as an effective catalyst to the engine cycle. It 
makes the energy transfer to the working substance faster. It is not called as true catalyst due to fundamental and practical 
constraints on its complete recovery by the time translation symmetry as well as rapid dephasing. Instead, it can be restored 
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before every interaction externally to imitate an effective catalytic use of coherence for enhancing energy transfer rate. SR 
enhanced thermalization is translated to work output of the photonic Otto engine, so that it is scaled quadratically with the 
number of atoms in the cluster. Such a scaling law is a quantum coherent effect, which cannot be realized with classical fuels, 
and hence exhibiting a profound difference of quantum fuels from their classical counterparts. 

Our proposed engine could work at higher temperatures, though still less than atomic and photon frequencies. In such 
a case our model should be generalized to Dicke Hamiltonian which includes counter rotating terms. Instead of preparing 
clusters at low temperatures, they could be at higher temperatures, including 7), > T c . Beyond a critical coupling strength and 
below a critical temperature, such a model predicts a superradiant phase. However, the critical coupling is in the ultrastrong 
interaction regime, g/(D > 1, and challenging experimentally. Instead, if we keep inducing coherence externally to the 
clusters, we find SR could still occur at a lower coupling strengths. Specifically, at g/(0 ~ 0.36, which is in SR regime, 
numerical simulations give 7 e fr,ss = T c + N 2 at finite T a ,T c < 1. The fluctuations around steady state values are less in the 
Dicke model case. Another advantage of inducing coherence is that the SR emission relies less on the critical time f crt , 4s 
which is required for the atoms to build up collective coherence inside the cavity. It is inversely proportional to the injection 
rate of the atoms in non-dissipative systems. By externally inducing coherence, we do not need to increase the injection rate 
to initiate SR earlier. Combination of noise induced coherence schemes with the coherence reservoir can be an attractive 
extension of superradiant QHE from the perspective of cost of coherence generation. 

In contrast to infinitely slow photonic Carnot engine, superradiant QHE in Otto cycle can produce finite power. It is 
necessary to make a dynamical simulation of the full cycle to examine the power output. We can predict the power scales 
linearly with the number of atoms in the cluster, as the thermalization time increases with the cluster population. 

We have recently examined a single multilevel atom as a quantum fuel for a photonic Carnot engine. The atom was 
assumed to have multiple ground states in coherent superposition, which are coupled to a common excited state by the cavity 
field. This effective model allows for analytical determination of the steady state photon number and related thermodynamical 
properties in Carnot cycle. It is found that the work output is scaled quadratically with the number of quantum coherent 
ground states. Rapid dephasing of the multilevel coherence as well as effective nature of the proposed level scheme and the 
vanishing power of the Carnot cycle make the benefits of quantum coherence in a single multilevel atom as a quantum fuel 
severely limited. An intriguing question is then if the quadratic enhancement in work harvesting capability, obtained both for 
a single multilevel atom as well as for many atoms as quantum fuels, is a fundamental limiting power law associated with the 
number of quantum resources; or if there can be more advantageous power laws, which can be translated into power output as 
well. Indeed, one can imagine combination of multilevel and SR enhancements. There are proposals that quantum entangled 
initial preparations for SR from multilevel systems could give faster than quadratic increase in intensity . Alternatively, SR 
emission in photonic crystals can also yield higher power laws. Our superradiant photonic QHE could be further improved 
by such extensions. 

Atomic clusters with effective catalytic quantum coherence can be used in other QHE cycles and systems. Their enhanced 
energy transfer rates make them fundamentally distinct quantum fuels from classical ones, and allow for technologically 
appealing power laws in the work output of QHEs. 


Methods 

Computational algoritm 

We run our simulations by using scientific python packages and some key libraries of QuTiP python. The numerical al¬ 
gorithm is prepared as one to one correspondence with the algorithm of the physical model. Inability to make adiabatic 
elimination of field mode ( g 2 N /k 2 1) and the requirement to take partial traces over atomic degree of freedom after each 
injection as well as the exponential increase of the Hilbert space dimension with the addition of atoms to the cluster set limits 
to the maximum number of atoms to be used to thermalize the working fluid. Again, the requirement to take partial traces 
over atomic degree of freedom makes the quantum Monte-Carlo trajectory methods unsuitable and force us to use full master 
equation approach. 

Moreover, the requirement of ordered injection of atomic clusters makes parallel programming over multi-cores almost 
impossible. On the other hand, we use multi-threading (parallel computing over a single core) to accelerate the code and 
reduce the run-time, though parallelising has no effect on the required dimension of the Hilbert space. 
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Effects of atomic decoherence 

The decoherence of the atomic clusters contributes to the dynamics as the spontaneous emission and the dephasing. The most 
general Lindblad dissipators reads: 

L{x)p = -^(2xpx' — x'xp — px'x), 

j N 

L(yi)p = ^£( 2 yiPyJ-yJyiP-pyJy,j ( 6 ) 

where x € {S . S z } if the contributions are collective and y ,■ S { cG , o‘} if the contributions are individual with a dissipation 
coefficient of y. Our simulations revealed that the superradiance emission profile prevails with slightly smaller effective 
temperatures for decoherence rates changing from 0 to 6k range and has the same behaviour as in in In Fig. 1(a) and 1(b). 
This result follows from the fact that the passage times of the atomic clusters are small for the atomic decoherence to kick in 
and the transfer rate of the quantum coherence is robust (~ N 2 ). The most detrimental effect is found to be pure dephasing 
and hence we consider pure dephasing rate as the main limiting factor to determine a most suitable experimental set up. 

In a typical microwave resonator with a frequency of co ~ 51 GHz, the atomic dephasing time 7? ~ 116 us leads to the 
dephasing rate y<f, ~ 10 4 Hz. The magnitude of the dissipation is about K ~ 10 3 Hz, therefore y,., ~ 1 Ok. Instead, if we 
consider an optical resonator system with a frequency of co ~ 350 THz, we have y</> ~ K ~ 10 8 Hz, though the coupling 
strength is about g ~ 220 MHz and therefore g ~ K. Thus, in microwave and optical resonator systems, the superradiance 
conditions are not satisfied. A circuit quantum electrodynamics set up, on the other hand, with co ~ 10 GHz, g ~ 100 MHz, 
y^/O) ~ 5 x 1 0 6 and k/co ~6x 1 0 4 seem to be most suitable modern resonator set up to satisfy the conditions for 
superradiant heat engine. 

Using circuit QED parameters, our numerical simulations give a range of 30 — 60 ns thermalization time for the number 
of atoms that we considered in the manuscript. The interaction time is taken to be ~ 0.1 ns. An optimization of the interaction 
time relative to atomic decoherence times require further examination to test if available control methods on atom-field cou¬ 
pling and uncoupling in circuit QED can be sufficiently fast. In addition, thermalization time is much larger than the typical 
time of adiabatic stages which should be larger than 1 /co ~ 0.1 ns. The power output of the engine is therefore limited by the 
thermalization time. 


Cost of coherence 

Let us consider a two level atom which is in its ground state |g). To create a coherent superposition state \\jf) = (1 /\/2){\g) + 
|e)), we may apply a square pulse with pulse duration T p and pulse amplitude E p . The area of the pulse is given by the relation 
A p = dEpZp/h with d being the magnitude of the dipole-matrix element such that the Rabi frequency is given by = dE p /h. 
We can approximate the pulse area with the tipping angle in the Bloch sphere, i.e, A ~ 6 which is necessarly equal to n/2 to 
obtain a coherent superposition state. If we denote the spontaneous emission rate of the two level atom with y, Fermi’s golden 
rule leads (n ~ 1) 

dl = lm^y 

CO 

where So is the vacuum permittivity, c is the speed of light and co is the resonant frequency of the field. We, then, obtain for 
the pulse amplitude that 


^ _ Tin I 
p 2x p Y 37T£o/!c 3 y’ 


and thus, the pulse intensity reads 

r _ C£ 0|p |2 _ 

b — -zr\ t p\ — 


2 1 Fl 24c 2 T 3 y 

For a beam with a width of 8 , the pulse energy is given by 

o k 2 co 2 8 2 

U p = I p n(8/2) 2 x p = hco- 


( 8 ) 


(9) 


( 10 ) 


96c 2 'Ey y' 

We make use of the relations co = 2nc/X and £ = kind, with A and C being the wavelength of field and the radial beam 
divergence, respectively. We obtain 


k- 11 n 2 

U p = hco— — — T ~ —h(o ~ 3hco, 

p 24 x P Y £ 2 3 


( 11 ) 
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where we set 1 /T ; ,y ~ 2 and £ ~ 0.5. For an N atom cluster, the cost of coherence is found to be f/ c ' ost =NU p . Thus, if we need 
m clusters to obtain superradiant themalization of the working fluid, the total cost of coherence is given by U c ost = ntU' cosX . 

If we use N = 2 atom ensembles m = 250 times to thermalize the working fluid, the maximum work output obtained from 
the superradiant heat engine is found to be W ou t ~ 0.35m ~ 10 25 J for a co/2n = 5 GHz resonator. This output is, then, at 
least three order of magnitude smaller than the cost J/ cos t = 1500m ~ 5 x 10 3 W ou t. 

It is crucial to distinguish the cost paid to generate the coherent atoms from the cost to maintain coherence in the working 
substance. The housekeeping cost to maintain coherence in working fluid can reduce the thermodynamic efficiency. In our 
case the atoms are the fuel and their coherence is the catalyst. Photon gas in the cavity is the working substance. Coherence 
does not need to be maintained in the working substance, photon gas, in the operation cycle of the engine. The generation 
costs associated with in the fuel are not included to the thermodynamical efficiency of engines. Accordingly we do not include 
the cost of coherence to the engine efficiency determination. 

Our calculation of cost of coherence produces a result larger than the work output of the engine and verifies thermody¬ 
namical constraints. One could consider a definition of round trip efficiency as another figure of merit for the relative energy 
harvested with respect to the total energy spent for the engine. From this perspective quantum coherence generation cost must 
be reduced to make the proposed engine more appealing for certain applications. 
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